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We have investigated whether the electron-phonon interaction can support a d-wave gap- 
anisotropy. On the basis of models derived from LDA calculations, as well as LDA linear-response 
calculations we argue that this is the case, for materials with buckled or dimpled Cu02 planes, for 
the so-called buckling modes, which involve out-of-plane movements of the plane oxygens. 
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The mechanism of high-temperature superconductiv- 
ity in hole-doped materials containing Cu02 planes re- 
mains a subject of vivid debate. □ A large amount of 
experimental data such- as superconductivity-induced 
phonon renormalizatipnsB, a large isotope effect away 
from optimal dopinga, r.and phonon-related features in 
the tunnelling spectraH show effects of the electron- 
phonon interaction. Density-functional (LDA) calcu- 
lations arc in agreement with much of this, but indi- 
cate that the electron-phonon interaction has insufficient 
strength (A s ~ lpu and leads to s-wave pairing. The 
experimental evidence that the symmetry of the paired 
state is dp with lobes in the direction of the Cu-0 bond, 
points to the Coulomb repulsion between two holes on 
the same copper site as the pairing agent. However, 
the scarcity of high-temperature superconducting ma- 
terials, as well as results of Hubbard- and t-J-model 
calculations!^, lead to the suspicion that something more 
than the Coulomb repulsion-ia-aeeded. 

It is therefore of interestJHHi3 to investigate whether 
the electron-phonon interaction, which gives a negative 
(attractive) pair-interaction, V ep (k, k') oc — \g (k, k') | , 
could support the observed singlet- and d-wave pairing 
with the gap-anisotropy: A (k) oc cos(ak x ) — cos (ak y ) 
sketched in Fig. [j]. For this to occur, and assuming a 
BCS-like zero-temperature gap-equation: 



2A(k') = -]>>(k',k; 



A(k) 



[e (k)-e^] 2 + A(k) 2 



the electron-phonon interaction must be large when A (k) 
and A (k') have the same sign, and small when they have 
opposite signs. The main pairing interaction is usually 
believed to be associated with the exchange of spin fluc- 
tuations, because it is remdsive and seems to peak for 
large momentum transferalJ. In the RPA, its form is 

where the band susceptibility xo (q,w) is supposed to 
be peaked near q = Ojl the basis of mod- 

els derived from LDA calculationsE3ii3, as well as LDA 
linear-response calculations^], we shall argue that this 
is the case for the so-called buckling mode, which for 
YBa 2 Cu 3 07 and q=0 is the 330 cm" 1 (hu -40 mcV) 
oxygen out-of-plane and out-of-phase mode. Although 
the interaction of electrons with this mode is proportional 



to the static dimple (YBa 2 Cu307; 5 — 7°) or buckle 
(doped La2Cu04 LTO phase), and hence expected to 
be small, Raman and neutron-scattering experiments, as 
well as LDA calculations, have shown it to be relatively 
strong. 




FIG. 1. Fermi surface of an optimally doped C11O2 su- 
perconductor (schematic). The energy bands are nearly flat 
around X and Y (slightly bifurcated saddle-points). For the 
experimentally observed d x i _ J/ 2-wave pairing, the supercon- 
ducting gap A (k) has opposite signs on the parts of the Fermi 
surface which are shown in full- and broken-line. The scattering 
of a Cooper pair at (k, — k) to one at (k', — k') is indicated. The 
pair-interaction V (k, k') which sustains this gap-anisotropy 
must be repulsive for momentum transfers q = k' — k connect- 
ing full and broken parts, and attractive for q's connecting full 
with full, and broken with broken parts. The first property is 
presumably provided by the Coulomb repulsion and the sec- 
ond, we argue, could be provided by the interaction between 
near-saddle-point electrons and buckling modes. 

Fig. 2 specifies our generic tight-binding model for 
the electronic band structure of a single Cu02-plane. In 
the upper part of the 2nd row are shown the four a- 
orbitals, \y) = 03 y , \d) = Cu x 2_ y 2, \x) = 02 x , and \s) = 
Cu s -(3 Z 2_i) , and in the lower part, the \d) orbital and two 
of the 7r-orbitals, \z) = 02 z and \xz) = Cu xz , seen from 
the side of the plane. When merely the pda and pdn hop- 
ping integrals t xd =t yd = i pd =1.6 eV and t z , xz =t z ^ yz =0.7 
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eV indicated in the 1st column are taken into account, the 
band structure and corresponding constant-energy con- 
tours, in Fig. 3, shown in the 1st columns arise: The 
cr-orbitals give rise to a bonding, a non-bonding, and 
an anti-bonding O a; -Cu a .2_ y 2-O y band plus a high-lying 
Cu s -(3 Z 2-i) level (full lines), and the 7r-orbitals give rise 
to two decoupled pairs of bonding anti-bonding bands 
which disperse in either the k x or k y direction (stippled 
lines). The orbital energies (in eV and with respect to 
the energy of the Cu a ,2_ J( 2 orbital) are indicated at the 
relevant points of the band structure. The conduction 



band is the anti-bonding pda band and the Fermi level 
will be close to the saddle-point at X (ir,0). In the 1st 
column, this saddle-point is well above the top of the 
pdir-band and is symmetric, in the sense that the abso- 
lute values of the band masses in the k x and k y directions 
are equal. This means that the constant-energy contour 
through the saddle-point is a square with corners at X 
and Y. This is shown in 1st row of Fig. 3. In the 2nd 
columns, the strong hoppings (t sp —2.3 eV) between the 
high-energy Cu s _(322_ 1 ) orbital and the 02 x and 03 y 
orbitals and the tiny 02 z and 03 z hopping are included. 
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FIG. 2. Eight-band Hamiltonian for a single C11O2 plane (2nd row) and synthesis of its band structure (1st row). 
The 1st column shows uncoupled pda (O a! -Cu :I ,2_ 3/ 2-Oy), spa (Cu s ) , and pdn (Cu xz -O z and Cuj, z -O z ) bands. In the 
2nd column, the coupling (t sp ) of the Cu s orbital to the pda band is included. In the 3rd column, also the coupling (t z d) 
between a- and 7r-bands induced by a static dimple (S = 7°) is included. Modulation of the latter gives rise to interaction 
with the buckling mode. All energies are in eV. 



This depresses the conduction band near (n, 0) and 
thereby increases the mass towards T (0, 0) and decreases 
it towards M (tt, tt) ; the saddle-point becomes asymmet- 
ric and the constant-energy contours now bulge towards 
r. With the flat part of the conduction band just strad- 
dling off the top of a 7r-band, even a weak dimple or 
buckle of the plane will introduce considerable hybridiza- 
tion between the a and pdn bands. This is seen in the 



3rd columns where we have turned on the weak Cn x 2_ y 2- 
O z hoppings (£ zd =0.24 eV oc sin(<5=7 )). Since the hy- 
bridization with the Cuzz-022 pdn band (but not with 
the lower- lying Cu yz -03 z band) vanishes at X (tt, 0), the 
saddle-point bifurcates away from X and towards T, once 
8 exceeds a critical value (~ 4° in YBa2Cu30y). The 
conduction band thereby becomes very flat in a region 
around X (and Y) extending towards T. 
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FIG. 3. Constant-energy contours in (k x , fc^-space (1st row), as well as in an (x,j/)-space where the constant-energy 
contours are hyperbolas (2nd row). The constant-energy contours shown are those of the second band from the top in Fig. 
2 and for energies close to that of the saddle-point at X (it, 0) , which is the Fermi level for optimal doping. The coordinate 
transformation is specified in the 4th column and the analytical expressions for the constant-energy contours are given in 
the 3rd and 4th rows. The functions d (e) , s (e) , and p (e) describe respectively the pder, spa, and air scatterings. In the 
first column s=oo, and p=0. In the second column s=0.6 and p=0. In the third column, s=0.6 and p=0.14. 



It is this near degeneracy with the top of a 7r-band of 
the saddle-point of the conduction-band, which has cr- 
character, which causes strong coupling to out-of-plane 
deformations. How near the a-ir degeneracy is, is first 
of all controlled by the energy of the Cu s _(3 Z 2_ 1 ) hy- 
brid and, hence, by the distance of apical oxygen from 
plane copper. Secondly, if the saddle-point of the a-band 
is suppressed below the top of the 7r-band and part of 
the Fermi-surface becomes 7r-like, then a static dimple 
or buckle will adjust itself through reduction of the ir- 
band width (oc t z ^ zx ) so that the 7r-band becomes essen- 
tially filled. Jjhis, we have learned from total-energy LDA 
calculationsEj. 

Although our tight-binding model has 8 orbitals per 
CuC>2 unit, the constant-energy contours, €j (k Xl k y ) = e, 
have simple analytical expressions. These are given in the 
bottom rows of Fig. 3 and are expressed in terms of func- 
tions, d (e) , s (e) , and p (e) , which describe the dpa, spa, 
and an scatterings. If we neglect the Cu xz and C\i yz or- 
bitals the constant-energy contours are merely hyperbo- 
las in a (x, y)-space defined by: x = h (1 — cos ak x ) and 
y = g (1 — cosaky) , and they are shown in the 2nd row 
of Fig. 3. The bottom rows give the explicit expressions. 
With this model we can now easily calculate 



Xi'faa;) = 27r/£&[/(e(k))-/(e(k + q))] 
x5(e (k + q) - e (k) - Twj) 

as a function of the band shape and the doping. For 
q ~ Q = (jr, 7r) , where the Coulomb repulsion is sup- 
posed to provide the pairing interaction, Xo vanishes 
when 6f is below the saddle-point and Hu> is less than 
~ fsaddic ~ £f (overdoping) . If cf is at or above a 
saddle-point and this is slightly bifurcated, then x'o 
has a large, broad peak around Q for hui larger than 
~ cf — e S addic- This peak is due to "nesting" of the 
flat band near X with the one near Y. This, we believe, 
gives rise to spin-fluctuations which provide the repul- 
sive, large-q part of the pairing interaction. For smaller 
q and near-optimal doping, x'o (li^- 1 ) with Tuo ^40 meV 
has high ridges which are caused by the "sliding" of a 
saddle-point along the Fermi surface. We shall now show 
that the bucling phonon takes advantage of this part of 
Xo and does not feel the peak at large q because here, its 

interaction \g (k, k')| 2 vanishes. 

The last point is illustrated in Fig. 4 where we have 
assumed a static dimple like in YBa2Cu307. The left- 
and right-hand sides show the electronic wave-functions 
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at X and Y, respectively. The Cu s _( 3z 2_ 1 ) character, ir- 
relevant for our qualitative argument, is neglected, but 
the weak Cu yz -03 z and Cua: Z -02 2 characters at respec- 
tively X and Y are included. As mentioned above, and 
as can be seen from the figure, there can be no Cu^- 
02 z character at X and no Cu yz -03 z character at Y. 
On the right-hand side, we have perturbed the wave- 
function at Y by the buckling mode with q = Q. This 
wave in the [1, Indirection increases/decreases the dim- 
ple of every other oxygen row and accordingly modulates 
the coupling t z d between each 02 z and its neighboring 
Cu^-^ orbitals. The electron-phonon matrix element 
('F (X)| 5 (Q) 'F (Y)} is seen to have contributions from 
hopping between orbitals, drawn with thick lines, only 
and it vanishes because S (Q) (Y)) is even and ('F (X)| 
is odd with respect to a y-axis through the 02's. 
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FIG. 4. Schematic illustration of why there is no elec- 
tron-phonon interaction for the bucling mode with Q = (71-, it). 

The strength of the electron-phonon interaction is 

d 2 q 



A = [irN (e F ) 



BZ 



BZA 



where 



N(e F ) 



dk 



|«(k)|- 



/FS BZA 
is the electronic density of states and 

2^/^{/[e(k)]-/[e(k + q)]} 



7(q,w) 



x6 [e (k + q) - e (k) - ?uj (q)] \g (k, k + q) f 
the phonon linewidth. For the cZ-wave channel with gap- 



anisotropy i (cos ak y — cos ak x ) 



-y, a factor (x — y) 



must be included in the integrand for N (e F ) , thus yield- 
ing Nd (cf) , and a factor (x — y) (x' — y') must be in- 
cluded in the integrand for 7 (q,w) . Hence for g con- 
stant, Xd would vanish. For the buckling mode our tight- 
binding model yields: 

3 (k,k') - 2t z 



(eF-e z )(ej?- 



y/(l - x) (1 - x') - y/(l - y) (1 - y') 



where for simplicity we have neglected the dispersion of 
the 7r bands. It is easy to see that this form vanishes for 
q = Q, for instance. To further illustrate this form we 




FIG. 5. (k — k')-dependence of the electron-phonon inter- 
action for the out-of phase buckling mode from the six-orbital 
model, k and k' were chosen on the YBa2Cu307-type Fermi 
surface, which lies between the two low hole doping con- 
stant-energy contours in the third column of Figure 3. 



show in Fig. 5 g (q) = 7 (q, 0) /xq (<lfi) , calculated for 
a band shape with a marginally bifurcated (extended) 
saddle-point slightly below ep. This g is seen to have 
the kind of dispersion required in order to support d- 
wave pairing, g (k, k') as given above, factorizes and the 
double-integral for may be expressed analytically and 
the result for uj — > is simply: A<j = c 2 o; _1 Nd (ef) , 
where c is the factor in front of the square parenthesis in 
the expression for g fk, k') . This result is positive definite 
for all dopings. For e F approaching a saddle-point at X, 
Nd — ► N, so that the strength of the electron-phonon in- 
teraction in this (i-channel is as large as in the s-channcl. 

Finally, in Fig. 6 we show the phonon spec- 
trum calculated ab initio with the linear-response LDA- 
LMTO method for the idealized infinite-layer compound 
CaCuC>2 doped with 0.24 holes per formula unitli3. In 
the calculations, this compound developed a static buckle 
(02 up and 03 down by 5=6°). The uppermost of the 
three phonon branches marked with spheres is the buck- 
ling mode and the numbers written on the modes are the 
mode A,j's. We see that the buckling mode dominates Xd 
and that it is large close to T along TX and vanishes at 
M (it, 7r). Numerically we found Xd = 0.3 and A s = 0.4. 
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FIG. 6. Phonon modes and the mode A<j calculated ab ini- 
tio with the linear-response LDA-LMTO method for CaCu02 
doped with 0.24 holes. 

We conclude that the linear electron-phonon interac- 
tion for buckled or dimpled CuCvplanes may support, 
but is hardly sufficient to cause, high-temperature super- 
conductivity based on d x i_ y i -pairing. The most impor- 
tant mode seems to be the buckling mode, because it 
modulates the saddle-points of the energy bands where 
the density of states is high and where the supercon- 
ducting gap is observed to be maximum. The electron- 
phonon interaction, for this mode, is small or vanishes 
for momentum transfers q near (7r,7r), and does there- 
fore not interfere with the <i-wave symmetry of the order 
parameter induced by spin-fluctuation exchange. 
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